$ 


rAO-A053  067  AIR  FORCE  GEOPHYSICS  '.AB  HANSCOM  AFB  MASS  F/§  17/9 

SHORT  ARC  REDUCTION  OF  RADAR  ALTIMETRY  COMPUTER  PROGRAM. (U) 

JAN  78  • HADGIGCORGC*  J TROTTER  NASA-P-57t70(8) 

UNCLASSIFIED  AF8L-TR-77-0170  NASA-€R-i81838  NL 


ADA053067 


Short  Arc  Reduction  Ot 


Radar  Altimetry  Computer  Program 


< 


tW\SA 

National  Aeronautics  and 
Space  Administration 


D D C 

pnOjDEI 

APR  24  1978 


SSMITE 


Wallops  Flight  Center 

Wallops  Island,  Virginia  23337 

OtSTFOVnON  STATEMENT  X 

AC  804  824-3411 

K jLppra?»d  far  public  icImm} 
Olsliljgofloa  UnUmitcd 

67% 


NASA  CONTRACTOR  REPORT  141434 


1 

i 


i 

' Short  Arc  Reduction  Of  Radar  Altimetry  Computer  Program 


George  Hadgigeorge 
and 

Jerry  Trotter 


Air  Force  Geophysics  Laboratory 
Air  Force  Systems  Command 
United  States  Air  Force 
Hanscom  AFB,  Massachusetts  01 731 


Prepared  under  Purchase  Order  No.  P-57270  (G) 


lUASA 

National  Aeronautics  and 
Space  Administration 

Wallops  Flight  Center 

Wallops  Island,  Vit^ginia  23337 
AC  804  824-3411 


KSTBlBUnbl^  statement  a 

^ Jtpproyed  for  publlo  ralecDMj 
pyuibatlOD  UnUmitad 


I 

i 


1 

i 

f 


TABLE  OF  CONTENTS 


DESCRIPTION 

PAGE 

BACKGROUND 

1 

INTRODUCTION  TO  SHORT  ARC 

APPLICATION  TO  REGIONAL 

GEOID  REDUCTIONS 

4 

INVESTIGATIVE  EXPERIMENTS 

CONDUCTED  TO  DETERMINE  THE 

FEASIBILITY  OF  THE  SHORT  ARC 

APPROACH 

8 

3.1 

Preliminary  Exercises  of 
the  Spheroidal  Multi quadric 

Model 

9 

3.2 

Computer  Simulation  of  the 
Recovery  of  Geoidal  Undula- 
tion Over  the  North  Atlantic 

14 

3.3 

Reductions  of  First  Available 
Satellite  Alti metric  Data  from 
Skylab 

20 

3.4 

Preliminary  Reduction  of 

GEOS- 3 Alti metric  Data  Over  the 
North  Atlantic 

27 

3.4.1  Unknowns  and  Constraints 

Used  in  the  Reduction 

28 

3.4.2  Results 

29 

MATHEMATICAL  SURFACE  MODELS 

35 

4.1 

Uses  of  the  Covariance 

Function  in  SARRA 

35 

4.2 

Spheroidal  Multiquadric 

Surface  Model 

37 

FORMATION  AND  SOLUTION  OF  NORMAL 
EQUATIONS 

39 

ANALYSIS  OF  RESIDUALS 

45 

iii 


SECTION  DESCRIPTION  PAGE 


7 

FINAL  RESULTS  OF  GEOID  SURFACES 
DETERMINED  BY  THE  SHORT  ARC 

APPROACH  UTILIZING  ALL  AVAILABLE 

DATA 

50 

7.1 

Unknowns  and  Constraints 

Used  in  the  Reductions 

50 

7.2 

Reduction  Results 

51 

7.3 

Residual  and  Profile 

Analysis 

56 

8 

PRE  AND  POST  PROCESSING  AUXILIARY 
PROGRAMS  FOR  SARRA  REDUCTIONS  OF 
SATELLITE  ALTIMETRY  DATA 

66 

8.1 

Pre-processor 

66 

8.2 

Post  SARRA  Plot  Programs 

70 

APPENDIX  A 

' 

73 

REFERENCES 

82 

iv 


1^  Aj* 


LIST  OF  FIGURES 


FIGURE  NO. 

DESCRIPTION 

PAGE 

1 

Illustrating  ground  tracks  of  every 
fifth  pass  of  approximately  200  passes 
of  GEOS-C  considered  in  hypothetical 
short  arc  determination  of  fine  struc- 
ture of  North  Atlantic  geoid. 

6 

2 

Map  showing  countours  of  Gravimetric 
Geoid  and  locations  of  selected  nodes 
employed  in  preliminary  testing  of  the 
Spheroidal  Multiquadric  Model. 

11 

3 

RMS  fit  - Grid  level  - Nodal  density. 

12 

4 

Contour  interval  - RMS  residual. 

13 

5 

Illustration  of  the  simulation  of  200 
satellite  passes  that  are  approximately 
over  the  North  Atlantic  Oceanic  region. 

16 

6 

Illustration  of  the  simulation  of  320 
satellite  passes  that  are  approximately 
over  the  North  Atlantic  Oceanic  region. 

17 

7 

Skylab  ground  tracks. 

21 

8 

Skylab  altimeter  data  pass  #4. 

23 

9 

Skylab  altimeter  data  pass  H - 
Puerto  Rico  Trench  only. 

24 

10 

Comparison  of  Gravimetric  Geoid  with 
Skylab  altimeter  data. 

25 

11 

Skylab  altimeter  data  pass  #9. 

26 

12 

Calibration  area  GEOS-3  ground  tracks. 

30 

13 

Geoid  contour  map  of  the  GEOS-3 
calibration  area. 

31 

14 

Comparison  of  GEOS-3  derived  Geoid  and 
NASA  Calibration  Area  Gravimetric  Geoid. 

33 

15 

Comparison  of  GEOS-3  derived  Geoid  and 
NASA  Calibration  Area  Gravimetric  Geoid. 

34 

V 


LIST  OF  FIGURES (coni^wued) 


FIGURE  NO. 

DESCRIPTION 

PAGE 

16 

Ground  track  of  satellite  passes 
over  the  North  Atlantic  Ocean. 

52 

17 

Geoid  contour  produced  by  the  reduction 
use  of  the  Covariance  Function  Surface 
Model . 

54 

18 

Geoid  contour  produced  by  the  reductions 
using  the  Multiquadric  Analysis  for  the 
Surface  Model . 

55 

19 

Residuals  from  pass  2531413. 

57 

20 

Residuals  from  pass  1222133. 

58 

21 

Residuals  from  pass  3040601. 

59 

22 

Geoid  profile  along  pass  2531413 
using  Multiquadric  Surface  Model. 

60 

23 

Geoid  profile  along  pass  2531413  using 
the  Covariance  Function  Surface  Model. 

61 

24 

Geoid  profile  along  pass  1222133 
using  the  Multiquadric  Surface  Model. 

62 

25 

Geoid  profile  along  pass  1222133  using 
the  Covariance  Function  Surface  Model. 

63 

26 

Geoid  profile  along  pass  3040601 
using  the  Multiquadric  Surface  Model. 

64 

27 

Geoid  profile  along  pass  3040601  using 
the  Covariance  Function  Surface  Model. 

65 

28 

North  Atlantic  data  rejection  limits. 

72 

Vl 


LIST  OF  TABLES 


TABLE  NO. 

DESCRIPTION 

PAGE 

1 

Basic  assumptions  underlying 

18 

various  simulations. 

2 

Key  results  of  simulations. 

19 

vii 


SECTION  1 

BACKGROUND 

In  1970,  AFGL  contracted  DBA  Systems  to  perform  a study 
of  various  reduction  approaches  to  satellite  altimetric  measurements 
of  the  oceanic  geoid  surface.  This  study  was  initiated  in  view  of 
the  then  upcoming  launch  of  the  GEOS-3  Satellite  which  was  designed 
to  carry  a radar  altimeter.  The  primary  objective  was  to  derive 
alternative  reduction  approaches  that  may  remove  the  necessity  for 
extremely  accurate  reference  orbits  of  GEOS-3  required  by  convent- 
ional approaches. 

Conventional  approaches  to  the  reduction  of  satellite 
altimeter  data  utilize  long  arc  orbit  integration  that  requires 
extensive  tracking  from  ground  based  trackers  to  maintain  extremely 
accurate  orbits.  A widely  held  premise  was  that  in  order  to  exploit  the 
satellite  altimeter  of  1 meter  accuracy  for  geoid  improvement,  it  was 
necessary  that  the  radial  component  of  satellite  position  be  known 
to  better  then  + 1 meter.  An  example  is  extracted  from  reference  11 
that  states  "Because  the  orbit  will  be  used  in  combination  with 
altimeter  measurements  in  the  data  reduction  process  leading  to  geoid 
improvement,  it  is  necessary  that  the  accuracy  in  S/C  height  as  cal- 
culated from  the  orbit  be  known  to  better  than  + 1 meter".  Problems 
arose  that  made  the  above  requirement  very  difficult.  Among  these 
problems  are  (a)  unreBolved  biaaea  in  a given  traaker  could  induce 


localized  systematic  errors  in  the  computed  orbit  which  could  be 
transferred  to  the  local  geoid,  (b)  unmodeled  perturbations  caused 
by  drag  radiation  pressure,  (c)  errors  in  long  arc  orbits  result 
from  integrated  effects  of  errors  in  the  adopted  geopotential  function 
and  (d)  most  import-ant  is  the  extensive  effort  and  cost  of  global 
tracking  networks  and  computer  processing  required  to  provide  such 
accurate  reference  orbits. 

In  view  of  the  above  difficulties,  our  investigation  was 
directed  toward  alternative  approaches  that  would  fully  exploit  the 
1 meter  altimeter  accuracy  for  geoid  determination  with  less  stringent 
requirements  for  orbital  accuracies.  This  led  to  the  investigation 
of  the  feasibility  of  utilizing  the  short  arc  technology  that  AFGL/ 

DBA  had  developed  in  previous  satellite  geodesy  programs.  Such  an 
approach  was  envisioned  to  involve  the  simultaneous  recovery  of 
state  vectors  defining  as  many  as  several  thousand  short  arcs  (sub- 
ject to  weak  or  prior  constraints)  along  with  the  mathematical  model 
used  to  define  the  geoid  surface.  The  estimated  orbit  accuracy 
requirement  was  + 100  meters,  which  is  easily  attainable  from 
routine  global  tracking. 

The  short  arc  approach  is  defined  as  orbital  arcs  no 
greater  than  one  fourth  revolution.  It  was  shown  in  Brom  {2967) 
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that  accuracies  of  integration  of  better  than  1 meter  can  be  attained 
in  short  arc  reductions  provided  that 

a)  spheriaal  harmonies  at  least  through 
(n^m)  = (4,4)  are  exeraised  in  the 
integration 

and 

b)  all  six  orbit-al  parameters 
(Xo,  Yo.  Zo,  Xo.  Yo,  Zo) 

at  mid-arc  are  free  to  adjust  to 
best  aaaomiodate  the  actual  orbit. 


SECTION  2 


INTRODUCTION  OF  THE  SHORT  ARC  APPLICATION 
TO  REGIONAL  6E0ID  REDUCTIONS 

The  initial  investigation  established  that  the  short  arc 
approach  was  feasible  and  offered  a flexible  means  of  utilizing 
altimeter  data  for  determining  the  geoid  surface  for  a wide  range 
of  applications.  The  initial  consideration  of  the  short  arc  method 
employed  spherical  harmonics  to  mathematically  define  the  geoid 
surface.  This  development  is  well  suited  to  global  representations 
of  the  geoid  undulations  and  has  been  implemented  into  a computer 
program  (sagg)  for  the  combination  of  Satellite  Altimetry  and 
Ground  Gravity  reduction.  However,  it  was  anticipated  that  the 
early  data  collection  phase  of  the  GEOS-3  program  would  concentrate 
on  altimeter  measurements  over  the  North  Atlantic  calibration  area. 
Additionally,  it  was  felt  that  investigations  of  fine  detail  over 
limited  localized  regions  would  require  an  extravagant  spherical 
harmonic  expression.  Therefore,  alternative  means  for  analytical 
representation  of  the  geoid  surface  were  explored.  The  most 
suitable  model  for  the  proposed  application  was  a derivation  of  the 
spheroidal  multiquadric  analysis  developed  by  Hardy  {1972).  In 
Hardy's  exercise  of  the  multiquadric  analysis,  the  nodes  {ooeffioient 
computation  point)  correspond  to  observed  data  points.  In  our  deriva- 
tion, the  nodes  do  not  necessarily  correspond  to  data  points,  are 
relatively  limited  in  number  and,  at  the  outset,  are  evenly  distributed. 


The  primary  objective  of  this  project  was  the  short  arc 
determination  of  the  North  Atlantic  geoid,  utilizing  the  spheroidal 
multiquadric  analysis.  This  led  to  the  development  of  the  computer 
program  SARRA  {Short-  Ara  Beduation  of  Radar  Attimetry).  It  was 
envisioned  that  the  altimeter  observation  would  cover  the  entire 
North  Atlantic  Ocean  at  spacings  of  approximately  a 1®  by  1®  grid 
with  an  initial  selection  of  nodal  points  at  5®  by  5®  grid  as 
illustrated  by  Figure  1. 

As  data  became  available,  it  was  evident  that  there  was 
a high  concentration  of  the  observations  over  the  calibration  area  with 
sparse  measurements  over  the  northeastern  portion  of  the  North 
Atlantic  Ocean  and  almost  none  in  the  southeastern  portion.  The 
results  of  the  reduction  indicated  a close  detailed  agreement  of 
the  dense  data  area  with  other  geoid  models  {Marsh,  Vincent,  Strange), 
but  the  geoid  in  the  sparse  areas  {northeastern)  seemed  tilted  toward 
the  more  dense  areas.  This  geoid  behavior  caused  concern  with  the  use  of 
the  multiquadric  analysis  when  the  data  is  non-uniform  and  sparse  in  cer- 
tain areas.  Consequently,  other  possible  surface  models  were  investigated. 

Upon  suggestions  by  Dr.  Donald  Eakhart  (AFGL).*  and  extensive 
theoretical  investigation  by  Dr.  George  Blaha  (DBA)**,  the  use  of  the 
covariance  function  was  introduced  as  a second  option  to  the  computer 
program  SARRA.  The  two  options  now  offer  flexibility  depending  on 
the  characteristics  of  the  available  data.  Both  options  are  demonstrated 

'*  Air  Force  Geophysics  Laboratory,  LGHanscotn  Air  Force  Base,  Ma. 

Personal  corrmmication. 

**  DBA  Systems,  Inc.,  Melbourne,  FL.  Personal  corrmmication. 
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FIGURE  1.  Illustrating  ground  tracks  of  every  fifth  pass  of 
approximately  200  passes  of  GEOS  C considered  in 
hypothetical  short  arc  determination  of  fine 
structure  of  North  Atlantic  geoid.  Indicated 
tracks  generate  approximate  5®x5°  cells.  Selected 
nodes  for  initial  representation  of  North  Atlantic 
geoid  by  means  of  spheroidal  multiquadric  functions 
are  located  nominally  at  alternate  corners  of 
5®x5®  cells  as  indicated  by  solid  dots.  Initial 
set  of  approximately  150  nodes  is  subsequently  to 
be  augmented  by  additional  nodes  at  locations 
indicated  by  residuals  from  initial  reductions. 


in  the  results  portion  of  this  report.  Details  of  the  multiquadric 
analysis  and  the  covariance  function  models  are  presented  in  Section 
4.1  of  this  report.  The  covariance  function  that  is  implemented  into 
SARRA  refers  to  a chosen  reference  ellipsoid. 

The  advantage  of  the  covariance  function  is  realized 
when  data  is  irregularly  distributed  with  some  surface  areas 
being  very  sparse  in  measurements.  The  covariance  function 
implies  some  apriori  knowledge  of  the  geoid  behavior  through 
the  spherical  harmonic  coefficients.  The  spheroidal  multiquadric 
analysis  is  completely  dependent  on  observation  data  and  reproduces 
a surface  model  as  the  best  least  squares  fit  of  the  observation 
data.  The  advantage  of  using  the  multiquadric  analysis  model  is 
realized  when  the  area  to  be  processed  is  covered  with  an  adequately 
dense  set  of  altimeter  measurements.  One  unique  feature  of  the  multi- 
quadric analysis  is  that  there  is  essentially  no  limit  on  the  size 
of  the  area  to  be  processed  and,  provided  enough  data,  one  can 
obtain  as  much  detail  as  the  data  provides  by  the  proper  selection 
of  nodal  points. 


SECTION  3 


P 


INVESTIGATIVE  EXPERIMENTS  CONDUCTED  TO  DETERMINE 
THE  FEASIBILITY  OF  THE  SHORT  ARC  APPROACH 


In  the  development  of  the  computer  program  SARRA,  a series 
of  investigations  were  conducted  for  the  purposes  of  determining 

a)  the  basic  feasibility  of  the  application  of 
the  short  arc  technology  to  geoid  recovery 
from  satellite  altimetry^ 

b)  the  most  suH-able  matti  models  for  represent- 
ing the  geoid  surface  in  regard  to  the  in- 
tended application^ 

c)  testing  the  overall  concept  with  simulation 
for  its  ultimate  accuracy  potentialy 

d)  formulation  of  the  overall  concept  into  a 
reduction  program  capable  of  routinely 
processing  the  observation  dato.  when 
GEOS- 3 became  operational 

and 

e)  develop  auxiliary  programs  for  the  pre  and 
post  processing  required  in  an  operational 
environment. 

The  key  steps  in  this  investigation  are  summarized  in  this  section. 
However,  the  full  detail  may  be  found  in  references  1 and  11. 

It  was  surmised  in  Brown  {1973)  that  the  spheroidal 
multi quadratic  model  would  be  well  suited  to  detailed  representa- 
tion of  the  geoidal  surface  over  such  local  regions  as  the  North 
Atlantic  (where  intensive  testing  of  GEOS-C  was  to  be  conducted). 
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In  such  applications,  the  model  has  the  virtue  of  relative  simplicity 
and  is  especially  attractive  in  the  flexibility  afforded  by  the  pro- 
cess referred  to  as  nodal  densifiaation.  This  is  an  adaptive  process 
wherein  the  original  set  of  regularly  spaced  nodes  {defining  the 
multiquadria  function)  are  supplemented  locally  as  needed  to  improve 
the  fit  over  irregular  areas  inadequately  modeled  at  the  outset. 

Despite  its  theoretical  attractiveness,  the  spheroidal 
multi quadric  model  had  to  be  implemented  and  tested.  The  next  section 
outlines  the  measures  taken  in  evaluating  the  model  under  a well 
controlled  experiment  with  data  typical  of  geoid  features. 

3.1  Preliminary  Exercises  of  the  Spheroidal  Multiquadric  Model 

In  order  to  ascertain  the  basic  adequacy  of  the  spheroidal 
multiquadric  model,  two  preliminary  numerical  tests  were  performed. 
The  first  test  was  to  determine  if  the  model  was  inherently  capable 
of  providing  a reasonably  good  fit  to  the  geoid;  the  second  was  to 
test  the  model  as  incorporated  in  SARRA  with  extensive  simulations. 
These  simulations  will  be  presented  in  the  next  section(Section  3.2). 

The  first  test  used  as  a data  base,  the  gravimetric  geoid 
of  a portion  of  the  North  Atlantic  produced  by  Marsh,  Strange  and 
Vincent  {i972).  A sample  of  about  250  spot  elevations  of  the 
geoid  were  extracted  from  the  contour  map  of  this  geoid.  These 
geoidal  heights  were  treated  as  if  they  were  direct  observations. 


The  multiquadric  model  was  then  fitted  by  least  squares  to  this  data 
set  using  an  initial  set  of  112  nodes  spaced  at  five  degree  intervals 
(see  Figure  2).  The  process  of  nodal  densification  was  invoked  in 
three  iterations  of  the  adjustment  to  introduce  fresh  nodes  in  re- 
gions leading  initially  to  residuals  in  excess  of  two  meters.  This 
led  ultimately  to  the  incorporation  of  22  additional  nodes,  several 
of  which  were  in  the  region  of  the  severe  undulation  over  the  Puerto 
Rico  Trench.  The  final  rms  error  of  the  fit  of  the  multiquadric 
model  turned  out  to  be  an  altogether  acceptable  0.75  m [the  precision 
of  the  digitization  of  the  contour  map  was  deemed  to  be  not  much 
better  than  0.5  m) . This  test  established  two  things: 

(a)  that  the  spheroidal  multiquadratic  function 
could  provide  an  accurate  representation  of 
regional  undulations  of  the  geoid, 

and 

(b)  that  the  process  of  nodal  densification 
provides  a valid,  effective  and  efficient 
means  for  extending  the  model  as  needed 
for  local  improvement  of  fit. 

Figure  3 is  the  contour  produced  by  the  reduction  and  Figure  4 
represents  contour  of  closeness  of  fit. 

With  the  fundamental  soundness  of  the  spheroidal  multi - 
quadric  model  firmly  established  as  a result  of  the  first  preliminary 
test,  steps  were  taken  to  incorporate  the  model  into  SARRA.  After 
the  revised  program  had  been  checked  out  on  DBA's  Xerox  Sigma  5 
computer,  a small  scale  simulation  was  performed  to  provide  a 
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FIGURE  2.  Map  Showing  Contours  of  Gravimetric  Geold  and  Locations 
of  Selected  Nodes  Employed  In  Preliminary  Testing  of 
the  Spheroidal  Multiquadric  Model 


Contour  Interval  = .5  meters 
RMS  Residual  = .75  meters 


tentative  indication  of  accuracies  to  be  expected  in  applying  the  model  to 
the  determination  of  geoidal  undulations  over  the  North  Atlantic. 

For  this  purpose,  a set  of  40  passes  was  generated,  half  of  which 
were  ascending  (SE-NW)  and  half  of  which  were  descending  (NE-SW). 

Nodes  (a  total  of  7i)  were  placed  at  alternate  points  on  a five  degree 
grid  and  a low  sampling  rate  of  one  point  ever  70  seconds  was  adopted 
(this  aovresponds  to  a spaaing  of  about  five  degrees  along  each  arc). 

A priori  one  sigma  constraints  of  5m  were  exercised  on  the  state  vectors 
for  each  of  the  40  arcs  and  a standard  deviation  of  Im  was  assigned  to 
the  altimetry.  Under  these  assumptions,  the  error  propagation  performed 
by  SARRA  indicated  that  rms  errors  in  the  recovered  geoid  could  be 
expected  to  range  from  a low  of  0.67  in  the  interior  of  the  region, 
to  a high  of  1.93m  near  the  boundary.  This  result  is  not  considered 
to  be  particularly  realistic  because  of  the  obviously  inadequate  num- 
ber of  nodes  exercised  in  the  model.  It  does,  on  the  other  hand, 
provide  some  indication  of  the  potential  power  of  the  method,  parti- 
cularly when  consideration  is  given  to  the  very  low  sampling  rate 
that  was  adopted. 

3.2  Computer  Simulation  of  the  Recovery  of  the  Geoid  Undulation 

Over  the  North  Atlantic 

The  successful  outcome  of  the  preliminary  test  made  it  logical 
to  proceed  with  large  scale  simulations  that  would  be  more  nearly  repre- 
sentative of  what  could  be  expected  from  the  GEOS-3  tests  to  be  conducted 

i 
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over  the  North  Atlantic.  For  this  purpose,  the  revised  version  of 
SARRA  was  converted  to  run  on  the  CDC-6600  computer  at  AFGL.  The 
simulations  performed,  consisted  of  two  runs  based  on  sets  of  320 
passes.  In  all  simulations,  the  nodes  were  spaced  at  regular  five 
degree  intervals,  the  total  number  being  126.  In  the  case  of  the 
200-pass  reductions,  the  altimeter  observations  were  assumed  to  be 
spaced  at  nominal  one  degree  intervals  along  each  arc  {this  oorres- 
ponds  to  a time  interval  of  14  seconds).  In  the  case  of  the  320- 
pass  reduction,  the  chosen  spacing  corresponds  to  0.625  degrees,  or 
a time  interval  of  8.5  seconds.  Schematic  layouts  of  the  ground 
tracks  for  the  two  sets  of  simulations  are  presented  in  Figures  5 
and  6.  Other  pertinent  data  concerning  the  distinguishing  assump- 
tions underlying  the  various  simulations  are  indicated  in  Table  1. 

In  Cases  3,  4,  5 and  6,  the  sigma  adopted  for  altimetry, 
namely  0.37m,  is  representative  of  what  would  be  expected  from  a 
data  compaction  process  based  on  the  following  considerations: 

(a)  an  original  sampling  rate  of  two  per  second, 

(b)  an  original  sigma  of  1.0m,  and 

(c)  use  of  a third  order  midpoint,  cubic 
polynomial  filter  employing  17  points. 

In  Cases  1 and  2,  no  data  compaction  scheme  was  considered  to  have 

been  exercised;  accordingly,  here  the  expected  sigma  of  the  raw 

altimeter  observations  was  exercised  at  the  aforementioned  data  rate 

of  one  point  every  14  seconds. 
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The  output  of  primary  interest  from  the  simulations  consists 
of  the  propagation  of  the  specified  errors  in  the  input  into  the  com- 
puted geoidal  heights.  Such  error  propagations  were  performed  for 
geoidal  points  corresponding  to  the  nadirs  of  the  data  points.  In 
due  course,  such  results  will  be  employed  to  generate  contour  maps 
defining  the  variation  in  the  error  of  the  recovered  geoid.  In  this 
section,  consideration  is  limited  to  the  results  indicated  in 
Table  2;  here,  the  high,  low  and  average  values  of  the  standard 
deviations  of  the  recovered  geoid  are  listed  for  each  of  the  simula- 
tions. 


TABLE  1.  Basic  Assumptions  Underlying  Various  Simulations 


CASE 

NUMBER 

OF  PASSES 

A PRIORI  ONE  SIGMA 
ORBITAL  CONSTRAINTS 
(Position/Velocity) 

ALTIMETRYa(m) 

1 

200 

20m/, 02  m sec”^ 

1.0 

2 

200 

15m/. 015  m sec  ^ 

1.0 

3 

320 

15m/. 015  m sec  ^ 

0.37 

4 

320 

15m/. 015  m sec  ^ 

0.37 

5 

320 

5m/. 005  m/sec  ^ 

0.37 

6 

320 

Im/.OOl  m/sec  ^ 

0.37 
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Table  2.  Key  Results  of  Simulations 


Case 

Number 
of  Passes 

One  Sigma 

Orbital  Constraints 
(Position/Velocity) 

Sigmas 

of  Recovered  Geo id 

KBB 

Low 

HESB^Bi 

1 

200 

20m/. 02  m sec"^ 

1.52 

1.34 

1.42 

■ 

200 

15m/. 01 5 m sec"^ 

.96 

.82 

.89 

B 

320 

20m/. 02  m sec"^ 

1.03 

.97 

1.00 

m 

320 

15m/. 01 5 m sec"^ 

.54 

.46 

.51 

320 

5m/. 005  m sec“^ 

.33 

.22 

.27 

1^1 

320 

Im/.OOlm  sec“^ 

.13 

.06 

.09 

The  results  indicate  that  general  accuracies  comfortably 
better  than  Im  {rms)  are  to  be  expected  from  the  reduction  of  320 
passes,  subject  to  a priori  orbital  constraints  on  the  order  of 
15m  (Case  4).  Because  orbital  accuracies  of  considerably  better 
than  5m  are  a reasonable  expectation  for  the  precise  reference  orbits 
ultimately  to  be  generated,  the  results  from  Cases  5 and  6 suggest  that 
geoidal  accuracies  on  the  order  of  a few  tenths  of  a meter  are  potent- 
ially attainable  through  the  application  of  the  Short  Arc  Method. 
Attainment  of  such  accuracies  in  practice  will,  of  course,  entail  the 
appropriate  application  of  nodal  densifi cation.  This,  in  turn,  will 
require  processing  of  additional  passes  to  maintain  a specified  level 
of  accuracy.  A reasonably  conservative  extrapolation,  in  our  view, 
is  that  with  the  exercise  of  some  50  to  75  well -placed,  additional 
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nodes  {raising  the  total  from  17b  to  200)  and  with  the  incorporation 
of  perhaps  an  extra  100  passes  {for  a total  of  over  400) , one  can 
expect  SARRA  to  produce  geoidal  accuracies  generally  better  than  0.5tn 
when  precise  reference  orbits  are  employed  in  the  adjustment  with  one 
sigma  a pr-*  jri  constraints  of  5m  {or  better). 

3.3  Reduction  of  the  First  Available  Satellite  Altimetry 

Data  from  Skylab 

The  first  satellite  altimetric  data  for  experimental  process- 
ing with  the  SARRA  computer  program  became  available  from  altimetric 
observations  made  by  Skylab  {SL-2).  Two  passes  were  processed  through 
the  SARRA  computer  program  for  the  purpose  of  testing  the  basic  data 
flow  procedures  and  to  evaluate  the  adequacy  of  the  program  in  reducing 
the  data.  Figure  7 shows  the  ground  track  of  the  two  passes  plotted  by 
the  DBA  CALCOMP  plotter.  Pass  9 begins  off  the  east  coast  of  the  United 
States  and  extends  approximately  15  to  20  degrees  below  the  equator. 

Pass  4 begins  off  the  east  coast  and  extends  over  Puerto  Rico. 

Mid-arc  state  vectors  were  computed  from  ephemeris  data 
provided  by  NASA  for  initial  orbital  estimates  of  each  arc.  Data 
editing  procedures  were  designed  and  performed  on  the  observation 
data  prior  to  SARRA  processing.  These  included,  primarily,  an 
automatic  rejection  of  gross  error  indicated  by  the  residual  altimetry 
parameter  in  the  SL-2  EREP  format  and  a three-sigma  criterion  edit 
after  a ninth  order  polynomial  data  fit. 
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The  approach  taken  in  this  experiment  was  to  process  the 
passes  through  SARRA  and  examine  the  behavior  of  the  minimized  resid- 
uals. Since  Pass  Four  extended  over  the  Puerto  Rico  Trench  vicinity, 
the  residuals  were  examined  particularly  for  undulations  in  that 
area.  The  residuals  after  SARRA  adjustment  are  shown  in  Figure  8. 

A sudden  drop  in  the  residuals  is  immediately  noticed  126  seconds 
from  epoch.  This  section  of  the  pass  was  extracted  and  processed 
again  with  only  the  observation  occurring  over  this  region.  The 
residuals  after  adjustment  are  shown  in  Figure  9.  This  residual 
profile  was  plotted  on  a figure  {Figure  lo)  reproduced  from  previous 
comparison  of  SL-2  versus  gravimetric  geoid  {Vincent,  Strange,  Marsh). 

The  residuals  of  Pass  Nine  are  presented  in  Figure  11. 

The  residuals  show  no  abrupt  changes  that  would  indicate  any  possible 
sharp  geoid  undulation.  There  is  a gradual  but  smooth  variation 
throughout  the  pass  of  approximately  7.0  meters  that  may  be  due  to  a 
combination  of  geoid  undulation  and  sea  state  variations. 

The  importance  of  this  experiment  was  simply  to  gain  exper- 
ience with  the  SARRA  program  with  real  observation  data  and  to  illus- 
trate the  sensitivity  of  the  program  in  detecting  geoid  undulations. 
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Result  from  DBA  SKYLAB  SL-2  altimeter  data  plotted  on  a reproduction  of  Figure  7 from  Marsh,  Vincent 
"DETAILED  GRAVIMETRIC  GEOID  FOR  THE  GEOS-C  ALTIMETER  CALIBRATION  AREA". 


3. 


Preliminary  Reduction  of  GEOS-3  Altimetry  Data  Over 
the  North  Atlantic  Ocean 


Some  preliminary  reductions  were  performed  with  GEOS-3 
altimetric  observations  concentrated  over  the  western  portion  of 
the  North  Atlantic  Ocean.  The  observational  data  was  provided  to 
AFGL  by  NASA  and  consisted  of  112  passes  {Figure  12).  These  data 
sets  provided  an  opportunity  to  test  the  basic  mathematical  models 
used  to  represent  the  geoid  surface.  The  results  of  these  reduc- 
tions were  presented  at  the  American  Geophysical  Union  {agu)  fall 
meeting  held  in  December  1976  at  San  Francisco. 

The  distribution  of  this  data  set  was  not  adequate  for  a 
geoid  undulation  determination  of  the  entire  North  Atlantic  and  was 
not  suitably  located  for  the  detailed  examination  of  such  fine  fea- 
tures as  the  Puerto  Rico  Trench.  However,  it  did  provide  adequate 
data  to  evaluate  the  performance  of  SARRA  computer  program  in  an 
operational  environment  and  demonstrate  the  potential  accuracy  of 
the  approach. 

The  short  arc  approach  to  these  reductions  may  be  visualized 
along  the  following  lines: 

(a)  reference  orbits  accurate  to  approximately 
20  meters  may  be  obtained  from  routine 
global  tracking; 

(b)  tlie  reference  orbits  are  divided  into  sub- 
arcs,  situated  over  oceanic  regions  and  are 
limited  in  length  to  no  more  than  a quarter 
revolution; 
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(c)  each  sub-arc  is  treated  as  an  independent  orbit 
with  the  epoch  at  mid-arc  having  a state  vector 
subject  to  a priori  constraints  consistent  with 
the  estimated  accuracy  of  the  reference  orbits: 

(d)  observational  equations  are  formed  from  satellite 
altimetric  measurements  for  each  sub-arc  for  the 
mathematical  model  chosen  to  represent  the  oceanic 
geoid; 

(e)  the  adjustment  simultaneously  recovers  the  coef- 
ficients of  the  geoid  surface  model  and  revised 
estimates  of  the  state  vectors  for  all  sub-arcs 
(wkich  oAe  unLimitzd  -Ln  numbeAJ. 


3.4.1  Unknowns  and  Constraints  Used  in  the  Reductions.  The  SARRA 
program  is  designed  to  determine  surface  coefficients  as  designated 
by  nodal  points  that  lie  within  the  boundaries  of  a dense  selection 
of  height  measurements.  In  addition  to  determining  the  coefficients 
of  each  nodal  point,  the  program  solves  the  six  orbital  parameters 
for  each  arc. 

The  orbit  parameters  were  assumed  subject  to  a priori 
constraints  of  20.0  meters  in  position  and  .005  m/sec  in  velocity. 

The  a priori  standard  error  was  assumed  to  be  1.0  meter  for  all 
altimetric  measurements. 

During  the  initial  preprocessing  of  the  GEOS-3  data, 
certain  characteristics  were  identified  for  editing  criteria.  The 
first  editing  level  automatically  examined  altimetry  measurements 
for  gross  errors  that  occurred  from  data  handling  procedures,  such 
as  tape  parity  errors  and  measurement  identification  problems.  The 


second  editing  procedure  was  to  examine  the  altimeter  measurements  for 
continuity  and  to  eliminate  abrupt  point  to  point  changes.  The  final 
editing  is  based  on  a three-sigma  criteria  when  compared  to  a polynomial 
smoothing  function. 

3.4.2  Results.  The  reduction  of  6E0S-3  altimetric  measurements 
over  the  North  Atlantic  included  data  from  112  passes.  The  surface 
grid  formed  by  these  passes  was  dense  in  the  western  portion,  and  very 
sparse  in  the  northern.  It  did  not  extend  to  the  southeastern  portion 
of  the  North  Atlantic.  Figure  22  illustrates  the  surface  geometry  of  the 
GEOS-3  ground  tracks.  The  figure  also  reflects  the  nodal  selection 
used  for  this  reduction.  Due  to  the  variations  in  data  densities 
from  one  area  to  the  other,  a proper  value  for  k {Equation  1,  Section 
3.2)  had  to  be  determined  for  representing  the  entire  North  Atlantic 
grid.  It  was  found  to  be  important  that  the  choice  of  6 = ka  bear 
a balanced  relationship  to  the  typical  spacing  of  the  nodes.  The 
more  closely  spaced  the  nodes,  the  smaller  the  logical  choice  for  6. 

This  allows  the  surface  in  the  vicinity  of  any  given  data  point  to 
be  determined  predominantly  by  those  kernel  functions  of  nearby  nodes 
and  yet  prevents  any  one  node  from  exerting  total  dominance. 

Figure  13  shows  contours  of  the  derived  geoid  undulations 
using  the  multiquadric  function  approach.  In  general,  the  entire 
geoid  agrees  favorably  with  the  Marsh  and  Chang  gravimetric  geoid, 
especially  in  the  western  portion  of  the  area  where  the  data  density 
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FIGURE  12.  Calibration  Area  GEOS-3  Ground  Tracks 


is  much  greater.  The  area  with  less  data,  such  as  the  northern 
part  of  this  region,  did  not  reflect  the  full  geoid  detail  shown 
by  the  gravimetric  geoid  used  for  this  comparison.  The  rms  of 
the  residuals  of  the  altimetric  measurements  after  the  adjustment 
was  1.8  meters. 

Geoidal  profiles  derived  from  two  passes  of  GEOS-3 
altimeter  data  have  been  compared  with  the  Marsh  detailed  gravi- 
metric geoid.  The  geoidal  comparisons  are  shown  in  Figures  14 
and  15.  As  can  be  noted  in  both  figures,  the  relative  shape 
agreement  for  both  passes  is  excellent. 


Area  Gravimetric  Geoid 


PASS  2160955 
NORTH  ATLANTIC  AREA 


SECTION  4 


MATHEMATICAL  SURFACE  MODELS 

As  stated  earlier,  the  model  originally  used  in  SARRA 
for  representing  the  oceanic  geoid,  was  a derivation  of  the  spheroidal 
multiquadric  analysis  developed  by  Hardy  {2972).  This  model  provides 
convenient  flexibility  when  processing  data  collected  over  extremely 
small  surface  areas  with  a dense  set  of  measurement  data.  It  is 
independent  of  spherical  harmonic  representation  of  the  geoid,  thereby, 
eliminating  the  requirement  to  solve  for  high  order  and  degree  spherical 
harmonic  coefficients.  However,  in  the  course  of  preliminary  data 
processing  of  GEOS-3  altimetric  data,  it  became  evident  that  the 
observation  data  was  very  sparse  in  certain  areas  of  the  North  Atlantic 
Ocean.  Therefore,  another  model,  namely  the  covariance  function  by 
Heiskamen  and  Moritz  was  investigated  and  implemented  into  the  SARRA 
computer  program  to  provide  more  flexibility  when  processing  sparse 
sets  of  data. 

4.1  The  Covariance  Function 

The  detail  derivation  of  the  covariance  function  as  applied 
in  the  SARRA  computer  program  may  be  found  in  Appendix  A of  this  re- 
port. The  expression  for  the  partial  derivative  matrix  is  defined 
here  for  the  formation  of  normal  equation. 
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The  covariance  function  D(ii/)  for  geoid  undulations  is 
obtained  by  averaging  the  product  N over  the  unit  sphere, 

= M{NN'};  (1) 

N and  in  this  expression  are  the  geoid  undulations  at  any  two 
points  separated  by  the  symbol  M indicating  the  average  over  the 
unit  sphere,  namely 

M {•}  = — qlff — (Odo, 


where  a represents  the  surface  of  this  sphere  and  do  is  the  element 
of  surface  area.  The  covariance  function  may  be  expressed  as  (Heiskanen 
and  Moritz,  1967): 


n=2 


P„  (cos  ip). 


Where 


dn 


R2  V 

= G"(n-1)" 


(2) 

(3) 


the  quantities  dp  are  called  "degree  variances"  for  geoid  undulation 

and  A?__  are  the  corrections  to  the  a priori  coefficients  of  the 
nm  — ^ 

spherical  harmonic  expansion  of  the  geopotential.  The  parameter  R 

is  the  mean  radius  of  the  earth  and  G is  the  mean  gravity  of  the 

earth's  surface.  By  i|<  = 0,  we  have  from  (2)  by  substitution 

00 

D(0)  = Z)  <l„  (2., 
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The  weighted  sum  of  covariance  functions  between  point 
and  a number  of  selected  nodal  points  Pj  gives  the  undulation 

y - 

where  cy  are  the  weights  which  must  be  determined  and  il^^yis  the 
spherical  distance  between  ?i  and  Py. 


4*2  The  Spheriodal  Multiquadric  Model 

The  expression  for  the  geoidal  model  expressed  in  terms 
of  the  spheroidal  multiquadric  functions  is  the  same  equation  (4) 
above.  The  only  difference  is  the  replacement  of  the  covariance 
function  D(i/j^y)  with  the  "kernal  function"  of  the  form 


ajjlXg  - Xyl  + ka]^  + ||Yg  - Yy|  + ka)^  * ||Z  . Z/l  + 


Where 

h-i  " ^s^  ■ h cos 
cos  (j.^  sin 


U) 
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And 


= Altimetry  measurement  of  the  geoid  point 

= Geocentric  radius  to  the  geoidal  point  Xg^,  Yg^,  Zg^ 

fea  = An  arbitrary  fraction  of  the  semi -major  axis  a 

Xy,  Yy,  Zy  = arbitrarily  specified  nodes. 

Xc;,  Yg;,  If-j  = geocentric  coordinates  of  the  ^th 
sub-satellite  point 
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SECTION  5 


FORMATION  AND  SOLUTION  OF  THE  NORMAL  EQUATIONS 


The  observation  equation  for  a single  pass  is  taken  from 
equation  (13)  of  Appendix  A and  written  as  follows: 


where 

= residual 

R|  = initial  estimate  of  the  radial  distance  to  the  sub-satellite  point 

d^  = small  quantity  to  correct  for  non-geocentric  direction  of  the 
altimetry  measurement 

dc^  = coefficient  matrix  of  nodes 

dx,...dz  = coefficient  matrix  of  orbit  parameters 


Here  the  subscripts  I and  / refer  to  the  sub-satellite 
point  and  the  node  of  a single  pass.  Later  the  subscript  k will 
be  introduced  to  refer  to  the  arc.  The  "kernal  function", 

is  selected  by  input  to  be  either  the  covariance  function 
or  the  multiquadric  function  to  represent  the  geoid  surface  as 
described  above  [.Section  4.1,  4.2). 

Equation  (1)  is  written  in  terms  of  the  matrix  of  partial 
derivatives  for  the  -i}^  altimetry  measurement  point  in  the  form. 


Where 


(2) 


Of,  = standard  error  of  the  altimeter  measurement 


Where 


taken  from  equation  (1)  as 


(3) 


.... 


n = number  of  nodal  points 


3(X,  ...Z)o 


14) 
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a 

I 
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The  parti als  with  respect  to  the  orbital  parameters  have 
been  refined  by  Bldha  {1976}  which  introduced  a faster  convergence 
in  the  adjustment  and  in  most  cases  removed  the  need  for  iterating 
the  solution.  The  detailed  description  of  the  parameters  in  (4) 
may  be  found  in  Appendix  A. 

The  discrepancy  term  is  defined  by: 

Et  = R.C  - n - + d-L  . 15) 

The  normal  equations  for  the  arc  are  written  as: 


where 


Nfe«fe=  Cfe 

- m 
Nfe  = ih  B, 


T nT 

= 1=1  B.C 


(5^  = correction  vector 


m = number  of  altimeter  measurements  in  the  fe  arc. 
The  normal  equations  generated  by  this  arc  may  be  written 


in  the  matrix  form 


\ "fe 

"fe  "fe 


Introducing  a priori  orbital  constraints  for  the  arc 


leads  to  the  following  system  of  normal  equations. 


i 

i 


Where 

= the  inverse  of  the  covariance  matrix  of  the 
a priori  values  of  the  state  vector. 


= the  difference  between  the  a priori  values  of 
the  state  vector  and  the  values  currently  being 
employed  as  current  approximations.  Initially, 
i may  be  considered  zero  for  the  first  iteration. 

The  expansion  of  equation  (8.)  is  now  introduced  for  the 
simultaneous  reductions  of  all  adjustable  parameters  for  all  arcs. 
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in  which 


(10) 


The  solution  of  the  normal  equations  (9)  involves  n unknowns 
(one  for  each  node)  plus  six  for  each  arc.  The  number  of  arcs,  s, 
in  a solution  may  grow  to  several  thousand,  thereby,  creating  a large 
system  with  rank,  r = n + 6s.  The  solution  to  this  large  system  of 
normal  equations  is  made  practical  by  virtue  of  exploiting  the  patterned 
characteristics.  Full  details  of  derivation  may  be  found  in  Brown  {1958) 
or  in  Brown,  Trotter  {1909).  The  computational  steps  follow. 

Compute  the  auxiliaries  for  the  pass 


Qfe  = (N  + ur\  n[ 

16, n)  (6,6)  (6,n) 

(n,n)  (n,6)  (6,n) 

Sjj  = N||j  Rfe 

(n,n)  ln,n)  (n,n) 

= c - Q 
(n,l)  (n,l)  in, 6) 


(11) 


(12) 


(13) 


16,1)  (6,6)  (6,1) 
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As  the  S|^  and  matrices  are  formed  for  each  pass,  they  are  summed 
into  the  matrices  S and  c by 


fe=; 


(15) 


and 


(16) 


fe=? 


The  solution  for  the  geoidal  coefficients  is  computed  by 
6-  S-1  c 


(17) 


(n,l)(h,n)(n,1) 

The  solution  6|^  for  the  adjusted  corrections  to  the  state  vector 
for  the  arc  is  computed  by 

^ h ) - Qfe  «fe-  (18) 

The  error  propagation  for  a given  geoid  point  may  be 
computed  by 

(19) 


where  the  Bp  partial  derivative  matrix  is  computed  for  the  corres- 
ponding 4),  X as  in  equation  (3). 
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SECTION  6 


ANALYSIS  OF  RESIDUALS 


Careful  pass-by  pass  examination  of  the  plotted  altimeter 
residuals  can  be  made  in  order  to  uncover  any  localized  systematic 
effects  signifying  insufficiently  modelled  geoid  undulations.  This 
is  especially  important  to  areas  corresponding  to  known  bathymetric 
and  surface  features  such  as  the  Puerto  Rico  Trench,  the  Mid-Atlantic 
Ridge,  sea  mounts,  etc.  Where  indicated  from  analysis  of  residuals, 
additional  nodes  can  be  introduced  at  appropriate  locales.  The 
process  of  nodal  densification  for  more  detailed  local  definition 

of  the  geoid  will  exercise  supplementary  observations  from  any 
additional  passes  that  may  be  available  over  that  area. 

Altimeter  residuals  will  reflect  not  only  unmodelled 
geoid  undulations,  as  just  discussed,  but  also  various  quasi- 
systematic  trends  attributable  to  ocean  dynamics  such  as  tides, 
wind  stress,  swells,  currents,  etc.  With  the  exception  of  currents, 
such  phenomena  are  ephemeral  or  cyclic  and  thus,  random  over  a 
sufficiently  large  number  of  passes.  Except  near  shore  lines,  their 
amplitudes  are  generally  less  than  one  meter  and  will  have  only  a 
slight  influence  on  the  geoid  recovery.  Although  the  residuals 
obtained  following  nodal  densification  contain  a wealth  of  infor- 
mation, it  is  beyond  the  scope  of  this  contract  to  subject  the 
residuals  to  such  a detailed  analysis. 
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This  section  will  present  the  procedures  that  may  be  pur- 
sued in  facilitating  the  ultimate  residual  analysis.  As  thousands 
of  passes  have  been  processed  with  the  computer  program  SARRA  over 
a course  of  months  or  even  years,  a very  stable  data  base  of  normal 
equations  and  coefficients  of  the  geoid  surface  will  be  established. 

As  new  observation  data  is  collected  {espeaially  in  regard  to  SEASAT  A), 
one  only  has  to  form  the  normals  of  the  new  passes  and  sum  them  into 
the  accumulated  normals  saved  in  the  data  base  in  order  to  update 
the  solution. 

Once  the  geoid  coefficients  become  sufficiently  stable 
(meaning  additional  measurement  makes  insignificant  changes  to  the 

geoid) i we  may  treat  the  geoid  coefficients  as  being  perfectly 
known  and  process  individual  passes  for  detail  residual  analysis. 

The  only  adjustable  parameters  are  the  six  state  vector  parameters 
(x^  7,  Z,  Xj  Y,  Z).  This  leads  to  a simple  adjustment  that  could 
be  processed  on  a mini -computer  or  a hardwired  microprocessor 
requiring  less  than  1000  words  of  memory.  It  could  even  be  envis- 
ioned as  a real  time  monitoring  processor  with  the  appropriate  com- 
munication and  transmission  interface. 

The  solution  of  the  single  pass  may  be  accomplished  as 

follows: 

the  observation  equation  is  written  as 


i 
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r . 

-c 


d (R  - r')i 
d (X,Y,Z,X,Y,Z) 


+ 


(1) 


d X 
d Y 
d Z 
d X 
d Y 
d Z 


[r  - r' + d - 


where  all  parameters  are  identical  to  those  described  in  Section  4. 
The  matrix  of  partial  derivations  of  the  altimeter  measurement 
is 


d (X,  Y,Z,X,  Y,Z) 


(2) 


the  discrepancy  term 


= 


r • + d • 


H-  + r. 

U. 


(3) 


The  geoid  parameter,  r^,  is  computed  from  the  geoid 
coefficients, 6,  equation  (18, Section  5)  saved  from  the  data  base  {previous 
SARRA  reduation)  and  partials,  , of  the  geoid  parameters  for  the 
measurement.  The  matrix  is  computed  and  formed  by  the  same 
process  as  in  Section  4.1,  equation  (2)  or  Section  4.2,  equation  (1), 
depending  on  which  model  was  used  to  represent  the  geoid  surface  in 
the  master  reduction.  The  solution  for  r^  is 

5 (4) 

(1,1)  (l.n)  (n,l) 
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where 


n = number  of  nodes  used  in  the  master  reduction 

and  the  partials,  B^,  are  computed  for  all  nodes  in 
the  master  solution. 

The  normal  equations  are  written  in  matrix  form  as 
N 6 = c (5) 

where 

m 

N = y bT  B . 

1=1 


c = 


6 = 


y bI  e, 

A.  A. 

1=1 


d X 
d Y 
d Z 
d X 
d Y 
d Z 


and 


m = number  of  observations  for  the  arc. 

In  order  to  make  the  above  solution  determinable  for  6,  the  orbital 
constraints  must  be  applied  as  in  Section  5. 

Equation  (5),  above,  is  modified  to  accommodate  orbital 
constraints  in  the  form 


(N  + w)  6 


= c - w E 


(6) 
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Since  the  value  for  C is  normally  zero  unless  some  correction  is 
known  in  regard  to  the  state  vector,  the  solution  for  6 becomes 


6 = (N  + wV  V 
(6,1)  (6,6)  (6,1) 


(7) 


By  substituting  S into  equation  (1),  we  solve  for  the  residuals, 
r^  for  all  points  on  the  arc.  These  could  be  output  to  plot  displays 
on  a CRT  type  monitor. 

The  above  process  requires  minimal  computer  core  requirements 

* • • • 

and  the  largest  matrix  required  is  (6,6)  for  the  N and  w matrix.  The 
processing  time  would  also  be  minimal. 
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SECTION  7 


FINAL  REDUCTIONS  OF  THE  NORTH  ATLANTIC  OCEANIC  GEOID 

The  major  objective  of  this  effort  was  to  develop  the 
appropriate  computer  software  for  reducing  altimetric  observations 
and  to  apply  this  software  in  the  determination  of  the  North  Atlantic 
oceanic  geoid.  The  initial  coverage  was  expected  to  be  similar 
to  that  depicted  in  Figure  1 and  become  more  densified  as  more 
observations  were  made  by  GEOS-3.  However,  the  availability  of 
data  did  not  provide  the  desired  coverage  over  the  entire  North 
Atlantic,  nor  the  density  in  key  areas  outside  the  calibration 
area.  This  data  set  provides  enough  coverage  (see  Figure  16)  for 
an  extensive  testing  of  the  capabilities  of  the  software  system  in 
detecting  and  modelling  detail  features. 

The  final  reductions  are  essentially  extensions  to  the 
preliminary  reductions  presented  in  Section  3.4 

7.1  Unknowns  and  Constraints  Used  in  the  Reductions 

The  same  set  of  nodal  points,  data  and  orbital  constraints 
were  used  for  both  the  spheroidal  multiquadric  and  the  covariance 
function  models.  The  number  of  nodes  used  was  47  and  the  number  of 
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passes  was  236.  This  resulted  to  a total  of  1463  unknowns  (ff  x 

passes  + nodes)  that  were  recovered  simultaneously.  The  constraints 

exercised  for  the  state  vectors  in  terms  of  north,  east  and  up  were 

a = 20.  0 meters 

n 

Og  = 20.0  meters 

= 20. 0 meters 

a*  = .010  meters/seaond 

n 

Og  = .010  meters/seaond 

Oy  = .010  meters/seaond. 

These  orbital  constraints  are  chosen  to  be  relaxed  enough  to  accom- 
modate expected  errors  in  the  generated  state  vectors.  The  alti- 
metric  observations  were  treated  as  subject  to  1.0  meter  error  {.one 
sigma)  for  the  adjustments  and  error  propagation. 

The  selection  of  nodal  positions  are  illustrated  by  Figure  16. 

7.2  Reduction  Results 

The  reduction  of  the  GEOS-3  altimetric  measurements 
include  data  from  236  passes  of  GEOS-3  over  the  North  Atlantic. 

The  surface  grid  formed  by  the  ground  track  of  these  passes  was 
very  dense  in  the  western  portion  and  sparse  elsewhere.  Figure  16 
illustrates  the  surface  geometry  of  the  grid. 
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It  has  been  experienced  (reference  11)  in  previous  appli- 
cations of  the  spheroidal  multiquadric  model  that  a proper  choice 
of  6 = ka  (Equation  1,  Section  3.2)  is  important  in  keeping  a bal- 
anced relationship  to  the  typical  spacing  of  the  nodes.  The  constant 
6 controls  the  degree  of  correlation  between  nodes  and  as  the  value 
of  6 decreases,  the  correlation  between  nodes  decreases  (see  refer- 
ence l).  The  present  reductions  varied  in  data  density  from  one 
portion  of  the  North  Atlantic  to  another,  thereby,  making  the  choice 
for  6 a compromise  between  the  dense  and  sparse  areas.  This  is  prob- 
ably the  primary  weakness  of  the  use  of  the  multiquadric  model.  How- 
ever, provided  uniform  grid,  adequate  density  and  the  proper  selection 
of  6,  the  geoid  surface  can  be  measured  to  a high  degree  of  local  detail. 

Figures  17  and  18  represent  the  geoid  contours  from  the  re- 
ductions utilizing  the  covariance  functions  and  multiquadric  models, 
respectively.  The  two  contours  generally  agree  to  about  one  to  two 
meters  in  the  very  dense  portion  of  the  North  Atlantic.  The  less 
dense  portion  of  the  North  Atlantic  differs  by  as  much  as  ten  meters. 

When  compared  to  the  Marsh  and  Chang  geometric  geoid,  it  appears  that 
the  reduction  using  the  covariance  function  agrees  more  favorably, 
although,  the  reduction  using  the  multiquadric  model  showed  better 
agreement  in  some  of  the  more  densely  grided  areas. 
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Contour  Produced  by  the  Reduction  use  of  the  Covariance  Function  Surface  Model 


The  geoid  contours  that  are  presented  in  Figures  17  and  18, 
are  automatically  generated  by  evaluating  the  nodal  coefficients 
obtained  from  the  adjustment  for  a set  of  coordinates  destributed 
throughout  the  observed  area. 

7.3  Residual  and  Profile  Analysis 

Once  the  coefficients  of  the  nodal  points  are  determined, 
the  residuals  are  computed  along  each  arc  for  evaluating  the  point 
to  point  fit  of  the  observations  to  the  adjustment.  As  stated  in 
Section  6,  these  residuals  reflect  certain  unmodelled  surface  fea- 
tures in  addition  to  the  random  noise  of  the  altimetric  observations. 
No  attempt  will  be  made  in  this  report  to  evaluate  the  meaning  of 
these  residuals  since  it  is  beyond  the  scope  of  this  contract. 

In  addition  to  computing  residuals,  the  SARRA  computer 
programs  computes  geoid  profiles  along  each  arc.  These  profiles  and 
the  residuals  are  saved  on  a magnetic  tape  or  disc  during  the  SARRA 
adjustment  and  on  option  may  be  plotted  for  a graphic  display.  The 
residuals  and  geoid  profiles  of  a few  selected  passes  are  presented 
in  Figures  19  through  27  as  examples. 
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FIGURE  21.  Residuals  from  Pass  3040601 


2531413 


o 


o 


CO  I— 

i-H  a; 

^ -o 

Jn  :? 

iO 

csj  a» 

u 

i/I  lO 
iA  M- 
(Q  U, 
O.  3 
(/) 
o> 

C U 
O •»“ 
f—  u 
fO  -O 
lO 
0)  3 

I—  O' 


o ^ 

t.  3 
O.  £ 

■o  o> 

C 

O 

0)  (/> 

O =3 


CM 

CM 


60 


1222133 


lO 


63 


FIGURE  25.  Geold  Profile  along  Pass  1222133  Using 
the  Covariance  Function  Surface  Model 
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Geold  Profile  along  Pass  3040601  Using 
the  Covariance  Function  Surface  Model 


SECTION  8 


PRE  AND  POST  PROCESSING  AUXILIARY  PROGRAMS  FOR 
SARRA  REDUCTION  OF  SATELLITE  ALTIMETRY  DATA 

Several  computer  programs  were  developed  under  this  contract 
to  perform  editing  {pre-proaessor)  prior  to  SARRA  reduction  and  dis- 
play {plots)  results  after  the  SARRA  reduction.  The  pre-processor 
served  to  edit  raw  altimetric  measurements  on  both  a gross  error 
and  on  a statistical  basis.  In  addition  to  editing  the  data,  the 
pre-processor  also  provides  a means  of  computing  a mid-arc  state 
vector,  given  the  sub-satellite  position.  This  was  required 
to  process  some  data  provided  by  NASA  where  only  (J),  X and  h were 
given.  The  plot  program  serves 

la)  to  plot  satellite  ground  tracks  of 
available  dat-a, 

(b)  plot  contour  map  of  resulting  geoid, 

{c)  geoid  profiles  along  a given  arc, 

and 

(d)  residuals  from  the  adjustment. 

8. 1 Pre-processor  and  Data  Noise  Filter 

Pre-processing  and  editing  {Prep)  of  GEOS-3  altimetric 
data  is  designed  as  a three-level  automatic  computer  editing  effort. 
The  raw  data  is  first  edited  for  "gross"  errors.  The  second  step 
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is  to  examine  the  data  for  continuity  and  eliminate  abrupt  changes. 

The  third  level  of  editing  rejects  data,  either  smoothed  or  unsmoothed, 
that  have  standard  deviations  greater  then  some  specified  criterion. 

The  editing  is  accomplished  in  three  steps, 

(1)  Qompaving  altimeti^ie  measurements  with  a 
reasonable  expected  range  of  possible  heights 
(h)  for  gross  error  such  as  parity  errors; 

(2)  continuity  is  detected  by  comparing  the  change 

in  the  measured  altimetry  (Lh)  between  successive 
data  points  with  the  corresponding  change  in  time 
(ht).  This  is  a measicre  of  the  altimetry  rate 
(H)  which  is  a smooth  function.  If  the  altimetry 
rate  exceeds  the  maximum  value y the  point  is 
rejected.  The  maximum  for  h may  be  estimated 
emperically  from  several  sets  of  GEOS-C  altimetry 
data. 

{3)  The  third  level  of  editing  rejects  data  based 
upon  an  input  sigma  criterion.  The  criterion 
will  be  based  upon  realistic  expected  accuracy 
of  the  GEOS-S  altimeter  measurements. 

The  criterion  chosen  in  step  3 should  be  a little  relaxed 
to  make  sure  that  points  are  not  rejected  that  truly  reflect  some 
surface  details  that  are  normal  sea  state  variations. 

In  addition  to  editing  measurement  data,  the  pre-processor 
can  compute  a mid-arc  state  vector  when  the  ephemeris  is  provided 
with  the  altimetric  data  tape  in  the  form  of  latitude,  longitude 
and  height.  The  following  computational  steps  are  implemented  in 
the  pre-processor. 

The  procedure  followed  in  computing  a mid-arc  state 
vector  for  each  pass,  was  to  compute  Cartesian  coordinates  from 
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the  provided  geographic  coordinates.  These  coordinates  were 
computed  by 


N * a/(l  - e^  sin2  <f>)  2 


’x' 

* (N+h)  cos  4>  cos  A 

Y 

= 

(N+h)  cos  sin  A 

.1. 

. [N(l-e2)  + hi  sin  0.. 

Where  a = the  earth  semi -major  axis 
e = eccentricity  of  earth 

A fifth  order  polynomial  fit  to  these  coordinates  was  performed 
in  order  to  reduce  the  effect  of  the  above  truncation  error  and 
to  provide  a means  for  computing  velocity  components. 

The  general  polynomial  expression  is  given  by 

■ *0  * *1*1  * *2‘i  • • • •„!? 

The  equation  written  in  matrix  form  for  all  i is  expressed  by: 


— “ 

1 tjtf  ... 

h 

1 t2t2  ... 

% 

• 

» 

♦ 

• 

• 

Xi 

1 t^tf  ... 

<■1 

• 
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Substituting  C for  the  vector,  B for  the  matrix  and  S for 
the  a^  vector,  yields 
C = BS  . 

The  least  squares  solution  [.with,  unit  weight  matrix)  for  the 
coefficient  vector  S is 

S = (B^  B)^  B^C  . 

The  order  polynomial  equation  is  evaluated  at  mid-arc  time, 
where 

*0  ^ - *j)  / 2 . 

The  velocity  component  at  X,  is  computed  from  the  first 
derivative  of  position  by 

This  procedure  is  repeated  for  the  Y,  Y and  Z,  Z 
components.  All  altimeter  observation  times  {ti)  are  initial- 
ized with  respect  to  mid-arc  by 

t.  = t.  - . 

The  velocity  components  corrected  for  the  earth  rotation  are 
computed  from 

w = *0  ■ ‘*'^0 

, y^  . Yq  + ?Xq 

io  * io  * 

to  obtain  inertial  components  (f  = earth  rotation  rate). 

\ 

\ 

\ 
\ 
\ 
\ 
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The  resulting  state  vector  representing  the  orbit  at 
mid-arc  (Xq,  Yq,  Zq,  x^,  and  s^)  is  input  to  the  SARRA  program 
with  the  appropriate  constraints. 

8.2  Post  SARRA  Plot  Programs 


The  plot  program  provides  the  capability  of  producing  satellite 
ground  track  plots,  geoid  contours  and  either  geoid  profiles  or 
residual  profiles.  Additionally,  the  land  features  superimposes  on 
the  ground  track  and  the  contours.  Initially,  the  land  outline  was 
read  from  maps  and  input  into  the  program  plots.  This  was  later 
replaced  at  AFGL  by  a magnetic  tape  containing  the  digitized  shore- 
line coordinates.  Each  plot  is  described  as  follows: 

(a)  G\ound  Tfuick  Vtot.  This  program  plots  the 
ground  track  of  satellite  passes  over  the 
North  Atlantia  from  a given  sub-satellite 
latitude  and  longtitude.  Additional  paral- 
lels and  meridians  may  be  selected  with 
program  input  parameters.  The  ground  track 
plots  may  be  contained  within  the  North 
Atlantic  data  boundaries  as  shown  in 
Figure  28. 

(b)  ContouA  PZoti.  Contours  of  either  the 
standard  deviations  or  the  actual  geoid 
heights  may  be  plotted.  The  standard 
deviations  are  computed  from  the  covari- 
ance matrix  from  the  SARRA  reductions. 

The  geoid  heights  are  computed  from  the 
surface  coefficients  saved  from  the 
SARRA  solution.  In  either  case,  the 
computations  are  performed  for  a grid 
of  surface  points.  The  generated 

grid  is  contained  within  the  bound-  ,■ 

aries  of  the  observed  geoid  and  the 

boundaries  of  the  North  Atlantic.  T 
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(c)  . The  altiTneter  observation 

residuals  (mejcu,ivizd  minuA  computed 
hetghtj  obtained  from  the  SARRA 
reductions  are  plotted  for  the  pur- 
pose of  visually  reviewing  measure- 
ment aharaateristias . The  residual 
profile  plot  provides  immediate  dis- 
play of  altimeter  measurements  noise 
and  such  oceanic  details  as  sea  state 
variations. 
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APPENDIX  A 


The  use  of  the  Covariance  Function  in  the  Short  Arc  Reduction 
of  Radar  Altimetry  Program 

The  following  derivation  of  the  covariance  function  and  its 
application  to  the  modelling  of  geoid  surface  in  the  SARRA  computer  pro 
gram  was  developed  by  Dr.  Georges  Blaha  of  DBA  Systems,  Inc.  under  a 
special  investigation.  The  results  of  Blaha's  investigation  is  re- 
produced as  an  appendix  for  detailed  reference  material  regarding 
Section  4 of  this  report. 


Covariance  Function  for  Geoid  Undulations 


The  covariance  function  for  geoid  undulations  (N)  is  derived 
from  the  basic  formula 


. n"  2 

where 


n 

o 


(1) 


R - is  the  mean  earth  radius  (R  = 6371  km), 

Y - is  the  mean  value  of  gravity  (y  = 979.8  gal), 

X - are  the  geocentric  latitude  and  longitude  of 
the  point  associated  with  N: 

the  other  overbars  indicate  that  we  are  dealing  with  "fully  normalized" 
harmonics.  We  have  (m  0): 


\ -'/i  n + 1 (sin  9), 

5^  cm  5)  - J2  (2„+i)  5), 
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(sin  ?)  and  (sin  9)  being  called  Legendre's  polynomials  and 
associated  Legendre  functions,  respectively.  Furthermore, 


where 

AC  «=AC  / -J  2 a +1  , (3s) 

no  no 


[ ^nm  I / 1 (n  + m) 

\ _ / V 2 (2  n + 1)*  (n-  m) 

I A S I 

I ™ I 


5 


(3b) 


and  where 


A C 


no 


is  the  correction  to  the  reference  C* 

no 

(to  be  mentioned)  in  order  to  obtain 


AC  * C . 
nm  nnr 


AS  = S ; 
nm  nm 


the  coefficients  C , C . 

no  nnr 

and  S's". 


S are  sometimes  called  "conventional  C's 
nm 


It  has  to  be  emphasized  that  equation  (1)'  is  valid  for  the  reference 
ellipsoid  having  the  same  mass  and  the  same  potential  as  the  geold.  This 
condition  is  fulfilled  by  the  "mean  earth  ellipsoid"  which  in  theory 
shares  with  the  actual  earth  two  additional  parameters  (cd  rotation  rate 
and  J2  ~Czo  )•  pertinent  constants  of  the  mean  earth  ellipsoid  in 

the  Geodetic  Reference  System  (GRS)  1967  are: 


C* 

so 

» - 1.0827  X 10“®, 

(4a) 

cio 

- + 2. 37  X lO'f 

(4b) 

Cto 

- - 6.1  X 10'®; 
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(4c) 

the  other  C*  are  essentially  zero  (already  C*  Is  very  small  and  is 
no  60 

sometimes  neglected). 

In  what  follows  we  shall  be  dealing  with  overbarred  coefficients 
related  to  their  conventional  counterparts  by  the  factors  shown  in  the 
equations  (3).  Adopting  the  mean  earth  ellipsoid  as  the  reference 
ellipsoid,  we  thus  have  in  agreement  with  the  equations  (4) : 


C|o  = -.484194  X 10 "3, 

Cjo  “ +.790  X 10"^ 

= -.17  X 10"®. 

60 

In  analogy  with  previous  notations 


A Cgo 

- C20 

0 

1 

A C^o 

= C 

40 

- C* 

AC 

= c 

- ^ . 

60 

60 

60 

(5a) 

(5b) 

(5c) 

(conventional  case),  we  write 

(6a) 

(6b) 

(6c) 


For  all  the  other  A we  have  essentially  (even  if  m = 0) 

AC  = C ; 
nm  nm 


all  A S are  in  fact  S themselves,  namely 
nm  nm  ' 


(7a) 


AS  “ S , 
mn  DID 


(7b) 


The  covariance  function  D (^)  for  geoid  undulations  is  obtained  by 
averaging  the  product  N N'  over  the  unit  sphere, 

D (♦)  - M (n  n'3  ; (8) 

M and  N'  in  this  expression  are  the  geoid  undulations  at  any  two  points 
separated  hy  the  spherical  distance  \|r  and  the  symbol  M indicates  the 
average  over  the  unit  sphere,  naaKly 
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a 

where  a represents  the  surface  of  this  sphere  and  do  Is  the  element  of 

surface  area.  It  can  be  shown  that  the  covariance  function  In  (8)  may 

be  expressed  from  (1)  as  follows: 

00 

D djj  Pn  (cos  )|f),  (9) 

n>=2 


where 


G*  (n-  1)2 


n 

Z(A  ^ + A S ^ ) ; 

nm  ma 

■ -0 


(10) 


the  quantities  d^  are  called  "degree  variances" 

If  \|f  «=  0,  we  have  from  (9)  : 

00 

■d<0).^ 

n-2 


for  geoid  undulations. 

(9') 


In  the  program  SARRA,  the  covariance  function  serves  to  express 
the  geoid  undulation  (Nj^)  at  a geoldal  point  as  a function  of 
parameters  cj  at  selected  nodes  P^: 


H 


1 


0 ('i'lp  * 

j 


(11) 


where  4^jls  the  spherical  distance  between  Pf  and  Pj.  The  overall 
adjustment  model  In  SARRA  can  be  written  as 

Hi  - Rf  - Tf  + dj. 
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where 


- Is  Che  radial  dlsCance  (from  the  coordlnace  origin  0)  Co 
the  satellite  point  S^, 

r^  - is  Che  radial  distance  to  the  subsatellite  point 

- Is  the  satellite  altimetry  observation  (distance  P^). 


The  small  unadjustable  quantity  d^  accounts  for  the  fact  that  0,  P^^,  and 
do  not  lie  on  a straight  line;  it  is  not  needed  explicitly.  The 
parameters  from  (11)  enter  into  the  adjustment  process  through  the 
model  equation 


ti  + 


where  rj^  Is  the  radial  distance  to  the  reference  ellipsoid  (it  is  computed 
from  (p^).  The  remaining  parameters  are  the  state  vector  components; 
these  parameters  are  always  weighted.  The  adjustment  model  in  SARRA  may 
thus  be  written  symbolically  as  follows : 

- R^  (state  vector  parameters)  - (c^  parameters)  - rj^  (9)  + d^.  (12) 

Vhen  performing  partial  differentiation  of  this  equation  with  respect  to 
the  parameters,  it  should  be  kept  in  mind  that  (p  in  rj^  (cp)  depends  to  a 
certain  extent  on  the  state  vector  parameters.  (See  reference  10  for  a 
detailed  discussion.) 


r^  = a/[l  +j~^2  sin^  being  the  usual  ellipsoidal  parameters 
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In  order  to  form  the  observation  equations,  from  (12)  we  deduce 


. 3 (R.  - r!) 

+ V.  = R°  - N?  - r!  + d + v ^ d (state  vec.  par.) 

1 t t 1 1 t o (state  vector  par.)  f ' 

0 

h Nj 


dTF-pr.)  ^ 

J o 


where  v^^  is  the  residual,  is  the  observed  altimetry  and  "o"  denotes 

the  approximate  (initial)  values  of  parameters  or  functions  of  such 

• • • 

parameters.  The  state  vector  parameters  are  X,  Y,  Z,  X,  Y,  Z in  the 
"Earth  Fixed"  (E. F. ) coordinate  system;  the  adjustment  in  SAREA  is 
effectuated  in  this  system.  From  (11)  we  have 


^ (c.  par.) 

J o 


~ I ('I' • • • # ® ('l'^  j) ^ • 1 » 


due  to  the  linearity  in  (11),  one  'can  take 


cj.o 


and  thus 


nJ  = 0, 


which  is  in  fact  used  in  SARRA. 


The  observation  eauation  is  then  written  as  follows: 

a 


dc. 


dX 


it 


(X,...Z)^  J 


+ (rJ  - \ + ‘*i)- 
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(13) 


With  self-explanatory  notations  this  equation  is  presented  as 


Vi  - a^  y + li. 


where 


(14a) 


(14b) 


(14c) 


The  term  li  is  sometimes  called  "discrepancy  term".  The  quantity  ti 
can  either  be  computed  as  in  (14c),  or  it  can  be  computed  as  the  distance 
OPi  from  the  coordinates  of  the  point  Pi  depicted  in  Figure  1.  This  Is 
In  fact  the  procedure  used  In  SARRA. 


, , , "ob8erved"geoid 


••.ellipsoid 


Figure  I 


By  Joining  together  Individual  observation  equations  along  one  short 
arc,  the  resulting  set  can  be  written  in  matrix  notations  as  follows: 


V - A )(  + L. 

Similar  procedure  would  apply  for  other  arcs  in  an  adjustment;  each  vector 
X would  now  contain  a different  set  of  state  vector  parameters  (besides 
the  cooiBon  set  of  Cj  parameters) . All  these  sets  would  finally  be  Joined 
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FIGURE  2; 


COVARIANCE  FUNCTION  FOR  GEOID  UNDULATIONS 
ASSOCIATED  WITH  DEGREES  TRUNCATION  N=2,  N=6. 
AND  N=:21 


into  one  large  system  of  observation  equations,  and  adjusted.  However, 
the  main  purpose  of  this  brief  exposition  of  the  SARRA  adjustment  model 
has  been  to  show  the  role  of  the  covariance  function  D (ijr) . We  have  seen 
how  to  compute  this  function  from  potential  coefficients.  We  shall  next 
discuss  some  practical  problems  associated  with  its  application. 

If  a set  of  potential  coefficients  (C's  and  S's)  is  given  complete 
through  a certain  degree  n = N,  the  covariance  function  D (f)  may  be 
constructed  for  varying  values  (between  0 and  Ztt)  according  to  the 
equations  (9*),  (9).  In  fact.  Figure  2 depicts  three  cases  of  the 
covariance  function  D (i)f)  ; they  correspond  to  the  truncations  N = 2, 

N = 6,  and  N = 21,  associated  with  a reasonable  set  of  coefficients 
complete  through  the  degree  and  order  (21,21).  The  figure  reveals  that 
the  greatest  contribution  to  the  value  of  the  covariance  function  comes 
from  the  lower  degree  and  order  coefficients;  there  is  a small  difference 
between  D (i^)  for  N = 6 and  N = 21  and  there  would  be  practically  no 
difference  if  N - 15  and  N ==  21  were  considered.  On  the  other  hand. 

It  appears  that  if,  for  example,  a given  (6,6)  model  were  to  replace 
the  reference  ellipsoid,  the  corresponding  covariance  function  of  new 
"undulations"  would  exhibit  sharp  and  distinct  features  which  would  vary 
substantially  for  different  degrees  of  truncation  (up  to  a relatively 
large  n = N).  It  Is  felt  that  such  an  approach  would  be  much  more 
sensitive  to  local  geoidal  features  Chan  is  the  use  of  Che  covariance 
function  D (i|r).  In  contrast  to  using  D (\|r),  such  a "modified  covariance 
function"  would  greatly  reduce  the  influence  of  distant  nodal  points  (e. g. 
when  f > 30*)  on  the  value  of  the  "undulation"  at  observation  points. 
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